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CN ■ Abstract 

c , 

The ideal neutrino detector at the neutrino factory should have a mass in the range of 
10 kton, provide particle identification to tag the flavor of the incoming neutrino, lepton 
qq ! charge measurement to select the incoming neutrino helicity, good energy resolution to 

reconstruct the incoming neutrino energy, and be isotropic to equally well reconstruct 
incoming neutrinos from different baselines (it might be more efficient to build various 
sources at different baselines, than various detectors). The detector should also be able 
to reconstruct neutrino event typically below 15 GeV. A detector with such quality is 
most adapted to fully study neutrino oscillations at the neutrino factory. In particular, 
^sO [ measurement of the leading muons and electrons charge is the only way to fully simul- 

taneously explore CP and T violation effects. We think that a magnetized liquid argon 
imaging detector stands today as the best choice of technique, that holds the highest 
promises to match the above mentioned detector requirements. We discuss also the op- 
timal neutrino factory energy and baseline between source and detector in order to best 
perform these studies. 
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The first generation long baseline (LBL) experiments — K2K §, MINOS §], OPERA (| 
and ICARUS j| || — will use artificial neutrino beams produced by the "traditional" 
meson-decay method, to search for a conclusive and unambiguous signature of the neu- 
trino flavor oscillation || observed in cosmic ray neutrinos Q. These experiments will 
provide the first precise measurements of the parameters governing the main muon dis- 
appearance mechanism. 

In contrast, a neutrino factory^ ^) is understood as a machine where low energy 
muons of a given charge are accelerated to high energy, and let decay into one electron 
and two neutrinos within a muon storage ring. 

The great physics potential of a neutrino factory comes from its ability to test in a 
very clean and high statistics environment all possible flavor oscillation transitions [[To], [ll| . 
This ability will provide very stringent information on all the elements of the neutrino 
mixing matrix and on the mass pattern of the neutrinos. 

In a 3 x 3-mixing scenario, the mixing matrix, which should be unitary, is determined 
by three angles and a complex phase. Neutrino factories will provide precise determina- 
tions of two angles and of the largest mass difference 1 10, PH . In addition, a test of the 
unitarity of the matrix could be performed |p"lf . 

a Based on an invited talk at the IX International Workshop on "Neutrino Telescopes", March 6-9, 2001, 
Venice (Italy). 
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Figure 1: View of the planned MINOS detector. The detector with a total mass of 5.4 kton 
should start data taking with the FNAL-NUMI beam in 2004. 

Apart from being able to measure very precisely all the magnitude of the elements 
of the mixing matrix, the more challenging and most interesting goal of the neutrino 
factory is the search for effects related to the complex phase of the mixing matrixfL^ . The 
complex phase will in general alter the neutrino flavor oscillation probabilities, and will 
most strikingly introduce a difference of transition probabilities between neutrinos and 
antineutrinos (CP-violation effects), and between time-reversed transitions (T- violation 
effects). Neutrino factories should provide the intense and well controlled beams needed 
to perform these studies. 

It should be stressed that the complete and comprehensive detection of T- and/or 
CP— violation effects is very difficult for terrestrial experiments, as it requires LjE v 
values simultaneously in the range of solar and atmospheric neutrinos. In addition, they 
require that the transitions v e ~ * v^, v e — > v^, — ► v e and — > v e be measured, a 
priori within the same experiment. It has been known that only in the case of the 
LMA solution to the solar neutrino problem, can one hope to look for effects related to 
the complex nature of the mixing matrix. We assume that this is the case. 

2 Detectors at the neutrino factory 

We briefly mention the kind of detectors that are currently envisaged in the context of 
the neutrino factory. 

2.1 Magnetized steel-scintillator sandwich 

This is the "traditional" neutrino detector, with a lot of experience gained with (though 
smaller) detectors like CCFR/NuTeV or CDHS. A detector based on iron has the ad- 
vantage of having a high density and to be easily magnetizable, for a "straight-forward" 
ji + 1 ' pT discrimination. It has sufficient granularity to only cleanly detect muons, and of- 
fers a rather poor discrimination and reconstruction of electrons and neutral-current-like 
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events (including hadronic decays of taus). The muon resolution is good and the jet en- 
ergy resolution is reasonable (typ. 80%/ \fEh). A minimum muon energy threshold (typ. 
4-6 GeV) is needed in order to separate the muon from other hadrons and the muon sep- 
aration from the jet is difficult to measure. The electron/hadron discrimination is rather 
poor. The angular resolution is determined by the transverse readout segmentation, which 
is in fact rather modest. The full volume of the detector has to be instrumented, and by 
nature of the sandwich, the readout is not very isotropic. Hence, a detector optimized to 
reconstruct "horizontal" events from an artificial neutrino beam coming from a well de- 
fined direction, will not at the same time provide good reconstruction of say atmospheric 
events coming from "above" and "below" . One considers as prototype for the neutrino 
factory the MINOS Q detector (see Figure [j]), that should reach a mass of 5.4 kton in 
2003. It is composed of 486 layers of 2.45cm Fe each, divided into two sections each 15 m 
long. The field has an average value of 1.3 T. Since detector construction is under way, 
we shall soon learn if this detector technology can be scaled to the larger masses currently 
envisaged for a neutrino factory, namely in the range of 40 kton. 

2.2 Large water-Cerenkov 

This is a well-proven technology (1MB, Kamiokande, SuperKamiokande) . The target ma- 
terial (Water) is cheap and only the surface of the detector (and not the volume unlike 
in the previous magnetized steel-scintillator sandwich) needs to be instrumented, hence 
this technology scales well to gigantic masses. The size will eventually be limited by 
water properties. A next generation 500 kton is under consideration |lj|. The reconstruc- 
tion of events is rather isotropic, and this kind of detector will certainly cover a broad 
physics program, including observation of atmospheric neutrinos, maybe solar neutrinos 
if energy threshold allows it, supernova neutrinos, and search for proton decays. The 
clear disadvantage of this kind of detector is that the reconstructed pattern is limited to 
"simple event topologies" , that are not really compatible with a detailed reconstruction 
of neutrino events with complicated final states. See Figure ||. The detection of the muon 
charge cannot be easily implemented in water and should rely on a downstream muon 
identifier. 




Figure 2: Simulated neutrino event from a 50 GeV storage ring. 
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2.3 Emulsion/target sandwich 

The ECC (Emulsion Cloud Chamber) technique currently envisag ed for OPERA § has 
been used successfully in the DONUT(l3) experiment, though in a very small scale com- 
pared to what is needed for a neutrino factory. The technique should be demonstrated at 
the 1 kton scale at the LNGS by the year 2005. Emulsions have a fantastic granularity 
(at the level of microns) and can be used to directly detect the kink in the decay of a 
charged tau lepton. In OPERA, emulsions will be used as very precise tracking devices 
to reconstruct pieces of tracks, whereas actual neutrino interactions will occur in a pas- 
sive target material (Pb). The non-alignment of the track segments reconstructed over a 
thickness of about 100 microns of two consecutive emulsion layers, will indicate potential 
decay kinks. A direct search for v e — + v T at a neutrino factory could be attempted if 
the charge of the tau can be detected to suppress the — > v T "background" . The most 
obvious disadvantages are the difficulty to scale this detector to very large masses, the 
potentially very large amount of scanning involved in a high-statistics experiment at a 
neutrino factory and possibly a severe background from charm decays produced in v e 
interactions. 
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Figure 3: Electronic liquid argon imaging of a cosmic ray induced shower. The overall drift 
time (horizontal axis) corresponds to about 40 cm of drift distance. The vertical co-ordinate 
is the wire numbering; around 40 cm are shown. 



2.4 Liquid Argon imaging TPC 

The liquid argon imaging technique provides fully bubble-chamber-like reconstruction 
of events (See Figure |3|). It is fully a homogeneous, continuous, precise tracking device 
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Figure 4: Side view during installation of one of the two cryostats that compose the ICARUS 
T600 detector. 



with high resolution dE/dx measurement and full sampling electromagnetic and hadronic 
calorimetry. Imaging provides excellent electron identification and electron/hadron sep- 
aration. Energy resolution is excellent (typ. 3%/y/E for e.m. showers) and the hadronic 
energy resolution of contained events is also excellent (typ. 30% /y/E~h). 

Like with the water Cerenkov detector, liquid argon detectors cover a broad physics 
program, including observation of atmospheric neutrinos, solar neutrinos, supernova neu- 
trinos, and search for proton decays. 

One disadvantage in the current implementation is the lack of magnetic field. One has 
in principle to rely on a down-stream muon spectrometer (that has however low threshold 
given the loss in Argon of 240 MeV/m). Magnetization was considered in the past |L5j, 
and a possibility to magnetize the large, multikton, volume of Argon is however under 
study |l6|]. This method is the only one, most easily scalable to multikton mass range, that 
would provide sufficient granularity to measure the charge of electrons (sec Section 5). 

Liquid argon imaging, though a priori more difficult to implement than say a mag- 
netized iron-scintillator sandwich, is becoming a mature technique, that has so far been 
demonstrated up to the 15 ton prototype scale. The next major milestone is the operation 
of the ICARUS 600 ton prototype (see Figure §). Its construction has been completed 
during 2000, in all its various components. The first technical run has started in March 
2001. Cosmic muon tracks have been seen in June 20010. The successful reaching of 
this milestone is very important, and after a series of other technical tests to be per- 
formed in the assembly hall within the summer 2001, the detector should be ready to be 



b See http://www.cern.ch/icarus/ for up-to-date information 
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transported to the LNGS tunnel by middle 2002. The physics program achievable with 
the T600 detector has been described in Ref.(5|. It covers the observation and study of 
atmospheric and solar neutrinos. 

3 The oscillation physics at the neutrino factory 

Neutrino sources from muon decays provide clear advantages over neutrino beams from 
pion decays. The exact neutrino helicity composition is a fundamental tool to study 
neutrino oscillations. It can be easily selected, since /i + — * e + v e i7^ and pT — > e~v & v ll can 
be separately obtained. 

At a neutrino factory, one could independently study the following flavor transitions: 

M~ — > e~ v e 

— > v e ^r e~ appearance (1) 

— * Vfj, disappearance, same sign muons (2) 

— » v T —>■ t~ appearance, high energy nu's (3) 

— > i7 e disappearance (4) 

— * Vfj, — > // + appearance, wrong sign muons (5) 

— > i7 T — ► r + appearance, high energy nu's (6) 

plus 6 other charge conjugate processes initiated from /i + decays. 

TTie ideal neutrino detector should be able to measure these 12 different processes as 
a function of the baseline L and of the neutrino energy E v ! 

Of particular interest are the charged current neutrino interactions, since they can 
in principle be used to tag the neutrino flavor and helicity, through the detection and 
identification of the final state charged lepton: 

vt.N — > £~ + hadrons DgN — » £ + + hadrons (7) 

We illustrate this in the case of a non-magnetized ICARUS-like detector. Figures ||, [| (?] 
and |§| show the reconstructed visible energy at the baseline L — 7400km normalized to 
10 20 /i's for each event class for a specific oscillation scenario with Am| 2 = 3.5 x 10~ 3 eV 2 , 
sin 2 023 = 0.5 and sin 2 2#i3 = 0.05. The different contributions including backgrounds for 
each event class have been evidenced in the plots. For example, in Figure |^, the different 
processes that contribute to the right-sign muon class are unoscillated muons, taus and 
background events. 

Hence, the ideal detector at the neutrino factory should possess the following charac- 
teristics: 

• Particle identification: the detector should be able to identify and measure the 
leading charged lepton of the interaction, in order to tag the incoming neutrino 
flavor. 

• Charge identification: the sign of the leading lepton charge should be measured, 
since it tags the helicity of the incoming neutrino. 

• Energy resolution: the incoming neutrino energy E v is reconstructed as E v = 
E( + Ehad, where Eg is the leading lepton energy and Ehad is the hadronic energy. 
Hence, detector with better energy resolution will reconstruct the parameter of the 
incoming neutrino better, and therefore the oscillation probability. 

• Low energy threshold: the reconstruction and identification should be fully ef- 
ficient for neutrino events below 15 GeV, as it is in this region that we expect the 
cleanest and most ambiguous signal from CP and T violation, as we will demon- 
strate in the Sections 9 to 12. 
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• Isotropic: one might want to perform various similar experiments at different base- 
lines. The probably most efficient way to achieve this is to build a large neutrino 
detector, isotropic in nature, capable of measuring equally well neutrinos from dif- 
ferent sources located at different baselines L. Because of the spherical shape of the 
Earth, sources located at different baselines L will reach the detector "from below" 
at different angles. 

It should be stressed that these features need to be implemented on detectors which 
have to be very massive. Indeed, the neutrino factory will clearly be more intense (by at 
least one order of magnitude) than current neutrino beams. However, the requirement 
to make precision measurements and the considered rather long baselines (neutrino flux 
scales as L~ 2 ), implies that detector in the range of 10 kt or more will be required. 

4 Detection of muons, electrons, taus and neutral cur- 
rents 

Clearly, the goal to identify and measure all possible final state leptons produced in 
charged current neutrino interactions as well as neutral current neutrino interactions 
imposes some constraints on the detector technology. 

The most stringent problem is related to the achievable granularity in a given detector 
configuration. Indeed, the requirement that detectors should have a fiducial mass beyond 
the kton range immediately brings in various choices and optimizations, since in first 
approximation, the total cost of the detector will depend on its mass times its granularity. 
It is therefore possible to opt for a poor granularity detector of very large mass, or a high 
granularity detector of a smaller mass. In between these two extremes one can have a 
continuous set of possible optimizations. 

When we consider the possibility to tag the outgoing lepton, we can subdivide the 
requirements as a function of the lepton type: 

• Muons: the detection of muons is straight-forward. In principle, one looks for 
penetrating particles. In practice, one has however to be careful of backgrounds 
(misidentification) coming from n^, decays and from charm scmilcptonic decays. 

• Electrons: the detection of electrons is harder. In principle, one looks for large and 
"short" energy deposition. In practice, one needs to carefully separate electrons from 
7T° conversions. Different levels of expertise have been developed in the field, yielding 
typically background rejections ranging from 1-2% (for the worst granularity) to 
better than 1CT 3 for the best granularity. 

• Taus: the detection of tau leptons is the hardest. One can attempt to look on 
an-event-by-event basis for the tau "kink" ; this however requires an extremely good 
reconstruction of the neutrino vertex. The required resolution has so far only been 
achieved with the help of photographic emulsions, which pose stringent constraints 
in scanning, and are strongly cost limited. Another possibility is to provide a 
"statistical separation" . In particular, 60% of the tau decays are into 1-prong, 3- 
prong or more hadrons and hence look like neutral current events (i.e. without final 
state leading electron or muon). 

As far as the detection of the lepton charge is concerned, here also the level of difficulty 
depends strongly on the type of lepton. In general, one has to rely on a magnetic analysis 
to measure the charge. 

• Muons: the measurement of the muon charge is relatively easy, since tracks produced 
by muons are long. One can cither envisage a fully magnetized target or rely on a 
down-stream spectrometer. 
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Figure 5: Visible energy spectrum for elec- 
tron events: v e CC (dashed line), v T and v T 
(dotted line) and v e CC (dot-dashed). The 
solid histogram shows the sum of all contri- 
butions. 



Figure 6: same as Figure || for right-sign 
muon sample: f„ CC (dashed line), v T and 
v T (dotted line) and meson decay background 
(dot-dashed). The solid histogram shows the 
sum of all contributions. 



30 CeV, L = 7400 km, id" \l* decays 




Wrong Sign fj, 

+ V T CC + background 
\CC 
V T CC 
Background 



Ami, = 3.5 x 1 CT ! eV J 
sin ! 0,, m 0.5 
sin 2 20,s = 0.05 




Figure 7: same as Figure ^ for wrong sign 

i / i \ .figures: same as b igure p for the no-lepton 

muon sample: CC (dashed line), v T and 



T> T (dotted line) and meson decay background 
(dot-dashed). The solid histogram shows the 
sum of all contributions. 



sample: v NC (dashed line), v T and v T (dot- 
ted line). The solid histogram shows the sum 
of all contributions. 
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• Electrons: the measurement of the electron charge is the hardest. One needs to 
measure significantly precisely the bending in the magnetic field before the start of 
the electromagnetic shower. Hence, the measurement is typically limited to a few 
X . 

• Taus: when the tau decays into a muon, then the measurement is easy. Otherwise, 
one needs to rely on a very high granularity and magnetized target, in order to 
identify the tau, its decay products and eventually reconstruct the mother charge 
from the charge of the decay products. 

5 Measurement of the electron charge in liquid argon 

We saw that liquid argon imaging provides very good tracking with dE/dx measurement, 
and excellent calorimetric performance for contained showers. This allows for a very 
precise determination of the energy of the particles in an event. This is particularly true 
for electron showers, which energy is very precisely measured. 

The possibility to complement these features with those provided by a magnetic field 
has been considered. Embedding the volume of argon into a magnetic field would not alter 
the imaging properties of the detector and the measurement of the bending of charged 
hadrons or penetrating muons would allow a precise determination of the momentum and 
a determination of their charge. 




Figure 9: Magnetized liquid argon TPC: simulation of the 2.5 GeV electron shower in liquid 
argon. The field has a strength B=1.5 T and is directed perpendicular to the sheet-plane. 

We have recently started studying the effect of the magnetic field on electrons (see 
Figure ^|). Unlike muons or hadrons, the early showering of electrons makes their charge 
identification difficult. We however found that the determination of the charge of electrons 
of energy in the range between 1 and 5 GeV is feasible with good purity, provided the field 
has a strength in the range of 1 Tesla. Preliminary estimates show that these electrons 
exhibit an average curvature sufficient to have electron charge discrimination better than 
1% with an efficiency of 20%. 

6 CP and T violation measurements 

In a three-family neutrino oscillation scenario, the flavor eigenstates v a {u = e, fi, r) are 
related to the mass eigenstates v[(i = 1, 2, 3) by the mixing matrix U 

v a = U ai pi (8) 
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(9) 



and it is customary to parameterize it as: 

C12C13 S12C13 s 13 e- zS 

U(9l 2 ,9 13 ,8 23 ,8) = ( -S12C23 - Ci 2 Sl3S23e Jl5 C12C23 - Si2Sl3S23e i<5 C13S23 
S12S23 - Ci 2 Sl3C23e I<5 -C12S23 - Si 2 Si 3 C23e l5 C13C23 

with sij — sin 9ij and = cos 6^ . 

We concentrate on transitions between electron and muon neutrinos. The oscillation 
probability is 

sin 2 A 23 s 2 2 Si3S23 + c 2 2 (sin 2 Ai3S 2 3 s 2 3 + sin 2 A 12 s 2 12 (l - (l + s 2 3 ) s 23 )) 
-^c 2 3 sin(20i 2 )si3 sin(26<23) cos<5 [cos2Ai 3 - cos2A 23 - 2cos(26<i2) sin 2 A i2 ] 

+ i C 2 3 sin(5sin(26'i 2 )si3sin(26'23) [sin2Ai 2 - sin 2A i3 + sin 2A 23 ] (10) 

where Ajk = Am^ k L / 'AE V (in natural units). This expression has been split in a first part 
independent from the phase 5, and in the two parts proportional respectively to cos S and 
sin S. To obtain the probabilities for — > v e and v e — > v^, we must replace S — ► —5, 
with the effect of changing sin S — ► — sin S and cos S — > cos 5. The term proportional to 
sin S is the CP- or T-violating term, while the cos S term equally modifies the probability 
for both CP-conjugate states. 

From this dependence, we see that a precise measurement of the v e — > oscilla- 
tion probability can yield information of the <5-phase provided that the other oscillation 
parameters in the expression are known sufficiently accurately. 

The dependence of the parameter 6 is a priori most "visible" in the energy-baseline 
range such that |Ai 2 | = |Ara 2 , 1 |Z/4E ! , ~ 1 and |A 23 | = \Aml 3 \L/AE v ~ 1. 

When |Ai 2 | <C 1 and |A 23 | ~ 1, we obtain 

P{v e ^v^) ~ ic 2 3{c 2 3A 2 2 + 2s 2 3 (sin 2 A 13 + sin 2 A23) (11) 



+2A 12 si 3 



sin(Ai 3 + A 23 ) cos 5 + (1 - cos(Ai 3 + A 23 )) sin S j 
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At even higher E v or smaller L, we further have, when both | A12 1 <C 1 and |Ai 3 |, | A23 1 <C 

P(v e -»i/ M ) =s ic 2 3{c 2 3A 2 2 + 2 s 2 3(A 2 3 + A 2 3 )+2A 12 s 1 3(A 1 3 + A23)cos ( 5}(12) 

and the dependence on the phase is only through cos 5. From this follows a degeneracy 
under the change of sign of <5. In this L and E v range, a precise determination of the 
oscillation probability can no longer determine the sign of 5. 

The behavior at various energies and baselines is explicitly shown in Figure [l^ for 
the two baselines L = 730 km and 2900 km as a function of neutrino energy E v . The 
probabilities are computed for three values of the <5-phase: 6 = (line), 6 = +tt/2 
(dashed), 5 = —ir/2 (dotted). The other oscillation parameters are: Am 2 2 = 3xl0~ 3 eV 2 , 
Amli = 1 x 10~ 4 eV 2 , sin 2 9 23 = 0.5, sin 2 6 12 = 0.5, and sin 2 26 13 = 0.05. 

A region corresponding to the oscillation of the "first maximum" is clearly visible on 
the curves. We define the energy of the "first maximum" as follows 

A 32 = - — ► E™ ax =^™™L 
32 2 v 2tt 

(1 x 1 97\ 

— > Er x (GeV) ~ Am 2 2 (eV 2 ) — L(km) (13) 
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Figure 10: Probability for v e — ► z/„ oscillations in vacuum for two baselines L = 730 km and 
2900 km as a function of neutrino energy E u . The probabilities are computed for three values 
of the 5-phase: 5 = (line), 5 = +ir/2 (dashed), 5 = —tt/2 (dotted). The other oscillation 
parameters are Am| 2 = 3 x 10~ 3 eV 2 , Am^ = 1 x 10~ 4 eV 2 , sin 2 9 2 3 = 0.5, sin 2 12 = 0.5, 
and sin 2 29 13 = 0.05. 



which yields E™ ax ~ 2 GeV at 730 km, E™ ax ~ 8 GeV at 2900 km and E™ ax ~ 20 GeV 
at 7400 km for Am| 2 = 3 x 10 -3 cV 2 . This energy corresponds to the point of maximum 
oscillation induced by Am§ 2 and coincides with the maximum when S = 0. It will be 
useful when we discuss the point of maximum sensitivity to the <5-phase. 

7 The rescaled probabilities 

In order to compare effects at different energies and various baselines, we define a "rescaled 
probability" parameter that allows a direct comparison of effects. Since to a good ap- 
proximation, the neutrino energy distribution at the neutrino factory behaves like Ef, and 
in addition, the neutrino flux scales like L~ 2 due to the beam divergence, we define the 
"rescaled probability" parameter p(i>a ^ vp;E v ,L) as 

E 2 

p{v a vp; E v , L) = P{v a ~* v - E v , L) x -| (14) 
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Figure 11: Rescaled probability (see text) for v e — > oscillations for a baseline L = 2900 km 
as a function of neutrino energy E v . The probabilities are computed for neutrinos in matter 
(full line) and in vacuum (dotted line), and for three values of the 5-phase: 5 = 0, 5 = +vr/2, 
S = — 7r/2. The other oscillation parameters are Am| 2 = 3 x 10 -3 eV 2 , Am|i = 1 X 10~ 4 eV 2 , 
sin 2 6» 23 = 0.5, sin 2 12 = 0.5, sin 2 29 n = 0.05. 



This behaviour is illustrated in Figure |TT| for a baseline L = 2900 km. We note that (1) 
it approximately correctly "weighs" the probability by the neutrino energy spectrum E„ 
of the neutrino factory spectrum; (2) it can be directly compared at different baselines, 
since it contains the L~ 2 attenuation of the neutrino flux with distance L; (3) p tends to 
a constant for E v — > inf , hence the high energy behavior can be easily studied. 



8 Propagation in matter 

Since neutrino factories will be associated to long baseline, it is not possible to avoid 
including effects associated to the neutrino propagation through the Earth matter. The 
simplest way to take into account these effects is to maintain the formalism developed for 
propagation in vacuum and to replace the mixing angles and the neutrino mass differences 
by "effective" values. 
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An important quantity for matter effects is D, defined as 

D(E v ,p) = 2V2G F n e E„ = 7.56 x l(T 5 eV 2 ( — P —){-^-) (15) 

gem a GeV 

where n e is the electron density and p the matter density. For antineutrinos, D is replaced 
by-D. 

There are two specific neutrino energies of interest when neutrinos propagate through 
matter: 

1. for D ;=y Am§ 2 , we reach for neutrinos the MSW resonance, in which the effective 
mixing angle sin 2 (29^) m 1. In terms of neutrino energy, this implies 

ET = COS2 ' 13Am ' 2 =g-cos2^ 13 Am 2 2 



2\/2G F n t 
1.32 x 10 4 cos 26<i3 Am^ 2 (eV 2 ) 
p(g/cm 3 ) 



in GeV (16) 



where £ res = (2\/2G_Fn e ) 1 . For density parameters p equal to 2.7, 3.2 and 
3.7 g/cm 3 one finds E r v es ~ 14.1, 12.3 and 10.7 GeV for Am| 2 = 3 x 1CT 3 eV 2 . 

2. for D > 2Am 2 2 , the effective mixing angle for neutrinos is always smaller than 
that in vacuum, i.e. sin 2 (26™\) < sin 2 (26> 13 ). In terms of neutrino energy, this is 
equivalent to 

E v > 2E r v es (17) 

3. these arguments are independent of the baseline L and depend only on the matter 
density p. 



9 Detecting the (!) phase at the NF 

In order to further study the dependence of the 5-phase, we consider the following three 
quantities which are good discriminators for a non-vanishing phase 5: 

1. A 5 = P(u e ^Vu,5 = +7T/2) - P(is e ^v„,5 = 0) 

The discriminant As can be used in an experiment where one is comparing the 
measured v e — > oscillation probability as a function of the neutrino energy E v 
compared to a "Monte-Carlo prediction" of the spectrum in absence of <5-phase. 

2. A CP (5) = P{v e -» 5) - P(D e - v„, 5) 

The discriminant Ac p can be used in an experiment by comparing the appearance 
of i/fj, (resp. Dp) in a beam of stored p + (resp. p~) decays as a function of the 
neutrino energy E v . 

3. A T (5) = P(v e -» 5) - P(y^ -» i/ e , 6) or A T (J) = P(P e -> ^, *) - -» ^e, 
The discriminant Ay can be used in an experiment by comparing the appearance 
of (resp. Dfi) and D e (resp. i/ e ) and in a beam of stored p + (resp. ^i~) decays as 
a function of the neutrino energy E v . 

Each of these discriminants have their advantages and disadvantages. 

The A<5-method consists in searching for distortions in the visible energy spectrum of 
events produced by the i5-phase. While this method can in principle provide excellent de- 
termination of the phase limited only by the statistics of accumulated events, in practice, 
systematic effects will have to be carefully kept under control in order to look for a small 
effect in a seen-data versus Monte-Carlo-expected comparison. In addition, the precise 
knowledge of the other oscillation parameters will be important, and as will be discussed 
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below, there is a risk of degeneracy between solutions and a possible strong correlation 
with the #13 angle at high energy. 

The Acp is quite straight-forward, since it involves comparing the appearance of 
so-called wrong-sign muons for two polarities of the stored muon beam. It really takes 
advantage from the fact that experimentally energetic muons are rather easy to detect 
and identify due to their penetrating nature, and with the help of a magnetic field, 
their change can be easily measured, in order to suppress the non-oscillated background 
from the beam. A non- vanishing Aqp should in principle be a direct proof for a non- 
vanishing (5-phase. This method suffers, however, from the inability to perform long- 
baseline experiment through vacuum. Indeed, matter effects will largely "spoil" Aqp 
since it involves both neutrinos and antincutrinos, which will oscillate very differently 
through matter. Hence, the Acp requires a good understanding of the effects related 
to matter. In addition, it involves measuring neutrinos and antineutrinos. The matter 
suppression of the antineutrinos will in practice determine the statistical accuracy with 
which the discriminant can be measured. 

Finally, the At is the theoretically cleanest method, since it docs not suffer from the 
problems of AS and Ac p. Indeed, a difference in oscillation probabilities between v e — > 
and — > v e would be a direct proof for a non-vanishing 5-phasc. In addition, matter 
affects both probabilities in a same way, since it involves only neutrinos. Unfortunately, 
it is experimentally very challenging to discriminate the electron charge produced in the 
events, needed in order to suppress the background from the beam. However, one can 
decide to measure only neutrinos, which are enhanced by matter effects, as opposed 
to antincutrinos in the Acp which were matter suppressed, and hence the statistical 
accuracy of the measurement will be determined by the efficiency to recognize the electron 
charge, rather by matter suppression. 



10 The LjE v scaling of the CP and T effects 

Regardless of their advantages and disadvantages, there is one thing in common between 
the three discriminants As, Acp and At: their behavior with respect to the neutrino 
energy E v and the baseline L. By explicit calculation, we find 



A s = --C13 sin 26»i2 si 3 sin 26*23 x 



(18) 



cos 2A13 — cos 2A23 — 2 cos 20i2 sin A12 + sin 2A12 — sin 2A13 + sin 2A23 
= — ^c 2 3 si 3 cos2Ai3 — cos2A 2 3 + sin2Ai2 — sin2Ai 3 + sin2A 2 3 
where for the second line we assumed for simplicity 6*12 = #23 = 7r/4, and similarly, 



C13S13 sin 26*12 sin 26*23 sin 8 



sin 2A12 - sin 2A i3 + sin 2A 



-C13S13 sin 26*12 sin 26*23 sin 8 



sin A12 sin A13 sin A23 



(19) 



As expected, both expressions vanish in the limit Am\ 2 — > where Am 2 3 — ► Am^. Also, 
as one reaches the higher energies, the terms Acp = At vanish as 



C13S13 sin 26*12 sin 2023 sin 5 Am 12 ttt sin Am 23 — — 



L 



4£„ 



L 



'4K 



— c i3 s i3 sin 20i2 sin 2023 sinJAm^Am^) 2 



L 

4K 



(20) 
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hence, in the very high energy limit at fixed baseline, the effects decrease as E~ 3 . That 
the effects disappear at high energy is expected, since in this regime, the "oscillations" of 
the various Ajk's wash out. 

The important point is that all expressions depend from some factors which contain 
the various mixing angles, multiplied by oscillatory terms which always vary like sine or 
cosine of A jk -terms (the terms in squared brackets in the expressions above). Hence, we 
expect the various discriminants to scale like Ajk oc L/E v . The sensitivity to the 5-phase 
will therefore follow the behavior of the oscillation probability, and we therefore argue 
that the maximum of the effect will occur around the "first maximum" of the oscillations, 
i.e. for E™ ax = Amj 2 L/2ir (see Eq. ©) g 

These considerations are strictly true only for propagation in vacuum. When neu- 
trinos propagate through matter, matter effects will alter these conclusions. We will 
however show that as long as the baseline is smaller than some distance such that the 
corresponding "first maximum" E™ ax lies below the MSW resonance neutrino energy 
E r J- s , the considerations related to the LjE v scaling are still largely valid. 

In what way does the matter effect alter the ability to look for effects related to the 
8 -phase? It is incorrect to believe that only the Acp discriminant will be affected by 
propagation through matter, since it is the only one to a priori mix neutrinos and an- 
tineutrinos. In reality, the "dangerous" effect of matter is to reduce the dependence of 
the probability on the <5-phase, and this for any kind of discriminant. 

In matter, we would for example write the At discriminant as 

| Ay | = (c^)V^ sin 26f 2 sin 20™ sin S rn [sin A™ sin A™ sin A™ 

- (c?3) 2 s?3 sm 26l i2 ^ 20 23 sin S [sin A™, sin A™, sin Agl (22) 

where because of our choice of Atom's, we have 8\\ ss 823 and S m w S. This implies 
that the <5-phase discriminants have a different structure that the terms that define the 
probability of the oscillation (i.e. the non (5-phase dependent terms). The discriminants 
are the products of sines and cosines of all mixing angles and of the Ajk's (see Eqs. ( |is| ) 
and (iL9|)). Because of this structure, their property in matter is different. 

The behavior for neutrino energies above the MSW resonance E r v es is determined by 
the fact that in this energy regime, 8 7 ^ } (E > E™*) — > it/2 and therefore — > and 
s i3 ~ ¥ 1- More explicitly, one can show that 



/ m\2 m 
\ C 13J °13 


oc E~ 2 




sin 2ff 2 \ 


~ sin 2023 = 


= const 


sin 20™, 


— > const. 




AMl x 


w AM 3 2 2 oc 


E v 


AMI, 


~ Aml 2 = 


const 



(23) 



c Strictly speaking, the maximum of the 5-phase sensitivity does not lie exactly at the "first maximum" as 
defined in Eq. From Eq. @, we expect the maximum to be "shifted" to higher values of L/E v , since 

it corresponds to the maximum of the term 

sinAi2sinAi3sinA 2 3 ~ Am^-r^rsin 2 ( Aml 3 — ^- ) (21) 



which has the functional form a; sin 2 x and, therefore, has its maximum shifted to higher values of x compared 
to sin 2 x. This small shift is smaller than the oscillation wavelength itself, and does not cause a major problem, 
since experimentally we will always have sufficient energy range to cover the full oscillation. 
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Therefore, 



A!p| cx E~ 2 sm5 



sin f Am 2 , -4- I sin 2 |AM?, ~ 



32 4£;,J olu v^" 13 4£, 



cx £" 3 (24) 



and one recovers a neutrino energy dependence identical to that in vacuum (see Eq. (BQ[)). 
Note also that the argument of the sine function AM^L/AE V is not small (i.e. the 
approximation sin a; ~ x is not valid). For our choice of oscillation parameters, the mass 
difference is approximately equal to AA/ 2 3 (eV 2 ) ~3x 10~ 4 x E^(GeV), and hence the 
dependence on the baseline is 



|A™| cx E~ 2 sm5 



sin 



fl,27Am 2 32 J^^ j sin 2 (3.8 x 10- 4 L(km))] (25) 



Hence, the discriminant will first be enhanced and then be suppressed by matter ef- 
fects. The maximum is found when the sine squared function reaches a maximum, or at 
approximately 4000 km under the assumption of high energy neutrinos. 

As anticipated, these discussions say that if one wants to study oscillations in the 
region of the "first maximum", one should not choose a too large baseline L, otherwise, 
matter effects will suppress the oscillation probability. This is even more so true, as it will 
be recalled below, that the magnitude of the effects related to the 5-phase are suppressed 
more rapidly than the oscillation. 

The simplest way to express the condition on the matter is to require that the energy 
of the "first maximum" be smaller than the MSW resonance energy: 

2V2G F n e E™ ax < Am 2 32 cos 29 13 (26) 

and, by inserting the definition of E r J lax = Am 2 2 L/2ir we get 



Em n x ^ 



'32 J 

7T COS 2#i3 7T COS 2#i3 



V2G F n e 2 x 1.27 x 7.56 x 10- 5 (eV 2 )p(g/cm 3 ) 
1.5 x 10 4 (km) 



p(g/cm 3 ) 



5000 km (27) 



To summarize, we find that the discriminants of the (5-phase all scale with 

LjE v . This is illustrated in Figure [l2[ where the rescaled At and Aqp discriminants 
are plotted as a function of the LjE v ratio. In the plots, sets of curves are shown for 
5 = +7r/2 and S — —ir/2. We see that the rescaled T-discriminant is (as expected) 
antisymmetric with respect to 6. For the shorter baselines (730 km and 2900 km) it is 
almost equivalent to the vacuum case (dashed curves). The 7400 km baseline yields a 
highly suppressed T-effect. The CP-discriminant has the same features, but is shifted 
with respect to zero due to matter enhancement. 

The most favorable choice of neutrino energy E v and baseline L is in the region of 
the "first maximum" given by (LjE v ) max ~ 400 for |Am§ 2 | = 3x 10~ 3 eV 2 . This leaves 
a great flexibility in the choice of the actual neutrino energy and the baseline, since only 
their ratio LjE v is determinant. 

Because of the rising neutrino cross-section with energy, it will be more favorable to go 
to higher energies if the neutrino fluence is constant. Keeping the L/E u ratio constant, 
this implies an optimization at longer baselines L. One will hence gain with the baseline 
L until we reach L max ps 5000 km beyond which matter effects will effectively reduce the 
dependence of the oscillation probabilities on the S-phase. 
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Figure 12: The rescaled At and Aqp discriminants (see text for definition) as a function 
of the L/E u ratio, computed for neutrinos propagating in matter at three different baselines 
L = 730 km, 2900 km and 7400 km, and also for propagation vacuum (independent of 
baseline). Three sets of curves are represented, corresponding to 5 = +n/2 and 5 = —ir/2 
(thick lines) and S = (thin lines). The other oscillation parameters are Am| 2 = 3x 10~ 3 eV 2 , 
Am^ = 1 x 10~ 4 eV 2 , sin 2 6» 23 = 0.5, sin 2 12 = 0.5 and sin 2 20 13 = 0.05. 

11 The correlation with #13 

We begin the discussion with the detection of the 5-phase with the method of the As 
discriminant. We recall that this method implies the comparison of the measured v e — > 
oscillation probability as a function of the neutrino energy E u compared to a "Monte- 
Carlo prediction" of the spectrum in absence of <5-phase. 

When searching for effects related to the 5-phase by comparing the measured v e — ► 

oscillation probability as a function of the neutrino energy E v to a "Monte-Carlo 
prediction" of the spectrum in absence of <5-phase, requires necessarily a precise knowledge 
of the other oscillation parameters entering in the oscillation probability expression. 

In particular, the knowledge of the angle #13 could be quite important. Indeed, the 
v e — * oscillation is primarily driven by the #13 angle and only thanks to a different 
energy dependence of the terms proportional to S than to those independent of S can one 
hope to determine #13 and 5 at the same time! 

This is however not true at high energy, when both | A.12 1 <§C 1 and | A.13 1 , IA23I <C 1. 
This can be explicitly shown for example for simplicity in the limit of small #13. The 
rescaled probability is in this case a constant: 

(Am 2 } 2 r 1 
p(i/ B -n/ M ) V 32 12j [l + 2Af.s 2 3 + 8jV Sl3 cosc?j (28) 

where M = ((Am 2 3 ) 2 + (Am^ 3 ) 2 )/(Am 2 2 ) 2 and N = (Am 2 3 + ATO^ 3 )/(Am 2 2 ). The 
absence of "oscillations" at high energy implies that a change of #13 can mimic a change 
of S. In practice, this implies that the best energy range to look for effects related to 
the (5-phase is close the "first oscillation maximum", i.e. E™ ax ~ 2 GeV at 730 km or 
E™ ax ~ 8 GeV at 2900 km. This implies that the detector should be able to reconstruct 
neutrino events at those energies with high efficiency, and low background. 



17 









= 7.5 GeV 




= 30 GeV 




Process 




L = 732 km 


L = 


= 2900 km 








10 21 \i~ 


2.5 x 10 20 n~ 




^ CC 




41690 




36050 


Non-oscillated 


v„ NC 




10700 




10300 


rates 


v e CC 




14520 




13835 




v e NC 




4266 




4975 


Oscillated 


Ve ~* Vji 


CC 


88 




50 


events (5 = vr/2) 


Vp, ~» V e 


CC 


258 




238 


Oscillated 


Ve ~* 


CC 


100 




54 


events (5 = 0) 


Vp V e 


CC 


385 




333 


Oscillated 


Ve Vfj, 


CC 


100 




55 


events (5 = —n/2) 


Vy, V e 


CC 


376 




330 



Table 1: Event rates for a 10 kton detector. The oscillation parameters are: Am 
3 x 10~ 3 eV 2 , Am 2 = 1 x 10~ 4 eV 2 , sin 2 (9 23 = 0.5, sin 2 la = 0.5 and sin 2 29 13 = 0.05. 



12 Two concrete examples at L=730 km and 2900 km 

In order to assess with concrete examples the use of the <5-phase discriminants, we consider 
the two baselines, with corresponding muon beam energy: 

• L=732 km, = 7.5 GeV, 10 21 muon decays 

• £=2900 km, E^ = 30 GeV, 2.5 x 10 20 muon decays 

Both examples were chosen to have the same L/E^. Because of the linear rise of the 
neutrino cross-section with E v , the factor 4 in muon energy between the 732 km and 
2900 km case, is "compensated" by an increase of intensity by the same factor in favor 
of the shorter baseline. 

The expected event rates are shown in Tables [l] and 0. 



12.1 Direct extraction of the oscillation probabilities 

From the visible energy distributions of the events, one can extract the oscillation prob- 
abilities. The visible energy of the events are plotted into histograms with 10 bins in 



energy. The v e 



oscillation probability in each energy bin i can be computed as 



Vi(v e 



Nj(w sfi) - N°(wsfi) 



m > P^*)iVP(e) 



(29) 



where Ni(wsfx) is the number of wrong-sign muon events in the i-th bin of energy, 
N®(ws(j,) are the background events in the i-th bin of energy, Ci(p^ > p™') is the ef- 
ficiency of the muon threshold cut in that bin, and Nf(e) is the number of electron 
events in the i-th bin of energy in absence of oscillations. The number of events corre- 
sponds to the statistics obtained from fi + decays. A similar quantity for antineutrinos 
Vi{y e — > v^) will be computed with events coming from pT decays. 

Similarly, the f M — > v e oscillation probability in each energy bin i can be computed as 



n Op 



N^wse) - N°(wse) 
e e (l - p C onf)N°{rsn) 



(30) 
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Table 2: Same as Table 1, but fi + decays. 



where Ni(wse) is the number of wrong-sign electron events in the i-th bin of energy, e e 
is the efficiency for charge discrimination, p CO nf the charge confusion, and N® (rs/x) is 
the number of right sign muon events in the i-th bin of energy in absence of oscillations. 
The number of events corresponds to the statistics obtained from fi" decays. A similar 
quantity for antineutrinos Viiy^ — > i7 e ) will be computed with events coming from /i + 
decays. 

These binned probabilities could be combined in an actual experiment in order to 
perform direct searches of the effects induced by the (5-phase. 

12.2 Direct search for T-asymmetry 

For measurements involving the discrimination of the electron charge, we limit ourselves 
to the lowest energy and baseline configuration (E^ — 7.5 GeV and L = 732 km), since 
we expect the discrimination of the electron charge to be practically possible only at these 
lowest energies. 

The binned Ay(i) discriminant for neutrinos is defined as 

A r (i) = Vtfa -> v e ) - Vi{v e -» !/„) (31) 

and a similar discriminant Ax(i) can be computed for antineutrinos. 

These quantities are plotted for neutrinos and antineutrinos for three values of the 
(5-phase (5 = +ir/2, (5 = and 5 = —n/2) in Figure [l3| The errors are statistical and 
correspond to a normalization of 10 21 muon decays and a baseline of L = 732 km. A 20% 
electron efficiency with a charge confusion probability of 0.1% has been assumed. The 
full curve corresponds to the theoretical probability difference. 

A nice feature of these measurements is the change of sign of the effect with respect of 
the change S — > — 5 and also with respect to the substitution of neutrinos by antineutrinos. 
These changes of sign are clearly visible and would provide a direct, model- independent, 
proof for T-violation in neutrino oscillations. 

In order to cross-check the matter behavior, one can also contemplate the CPT- 
discriminants defined as 

A C pt(«) = Vi{vn -> »e) - Vi(y e -> Vfj) 

K CPT {i) = Vi{v e ->!/„)- Vifa -> u e ) (32) 
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These quantities are independent from the ^-phase and probe only the matter effects. The 
change of sign of the effect with respect to the substitution of neutrinos by antineutrinos 
is clearly visible. 

It should be however noted that in the case of the CPT discriminant, the statistical 
power is rather low, since this measurement combines the appearance of electrons (driven 
by the efficiency for detecting the electron charge) and involves antineutrinos, which are 
suppressed by matter effects. Hence, the statistical power is reduced compared to the 
T-discriminant . 



12.3 Direct search for CP-asymmetry 

In the direct search for the CP-asymmetry, we rely only on the appearance of wrong-sign 
muons. We compare in this case the two energy and baselines options. 

The binned Acp(i) discriminant for the shortest baseline L — 732 km, — 7.5 GeV 
and longest baseline L — 2900 km, £L = 30 GeV (lower plots) for three values of the 
(S-phase (8 = +tt/2, 8 = and 8 = — tt/2) are shown in Figure [l3|. The errors are 
statistical and correspond to a normalization of 10 21 (2.5 x 10 20 ) for L = 732(2900) km. 
The full curve corresponds to the theoretical probability difference. The dotted curve is 
the theoretical curve for 8 = and represents the effect of propagation in matter. 

As was already pointed out, the Acp does not vanish even in the case 6 = 0, since 
matter traversal introduces an asymmetry. At the shortest baseline (L = 732 km), these 
effects are rather small. This has the advantage that the observed asymmetry would be 
positive for 8 > 0, but would still change sign in the case 8 ~ —tt/2. In the fortunate case 
in which Nature has chosen such a value for the 5-phase, the observation of the negative 
asymmetry would be a striking sign for CP-violation, since matter could never produce 
such an effect. 

For other values of the <5-phase, the effect is positive. It is also always positive at the 
largest baseline L = 2900 km, since at those distances the effect induced by the <5-phase 
is smaller than the asymmetry introduced by the matter. 



12.4 Comparison of two methods 

The binned Ap(z) and Acp(i) discriminant can be used to calculate the \ 2 significance 
of the effect, given the statistical error on each bin. We compute the following x 2 's: 

a = v (A T M)-A T M = 0)) 2 (A t M)-A t M = 0)) 2 
XT Y ^(AtM)) 2 <t(AtM)) 2 1 ' 

where a(Ax(i, 8)) is the statistical error in the bin. Similarly, the x 2 of the CP-asymmetry 
is 

2 _^ (A C pM)-A cp M = 0)) 2 

Y o{Acp{i,o)y 

We study the significance of the effect as a function of the solar mass difference Am^, 
since the effect associated to the (5-phase will decrease with decreasing Ato 2 ^ values. 
We consider the range compatible with solar neutrino experiments, 10~ 5 < Am^ < 
10- 4 eV 2 . 

The exclusion regions obtained at the 90%C.L. (defined as A% 2 = +1.96) in the 8- 
phase vs Am\ x plane are shown in Figure [l4| A 20% electron efficiency with a charge 
confusion probability of 0.1% has been assumed. The normalizations assumed are 10 21 
and 5 x 10 21 muon decays with energy = 7.5 GeV and a baseline of L = 732 km. 



20 




Figure 13: Direct T- and CP-violation: Binned At (i) (left) and Acp(i) (right) discriminants 
for neutrinos (upper plots) and antineutrinos (lower plots) for three values of the <5-phase: 5 = 
+7r/2, 5 = and S = —tt/2. The errors are statistical and correspond to the normalizations 
given in the text. A 20% electron efficiency with a charge confusion probability of 0.1% has 
been assumed. The full curve corresponds to the theoretical probability difference. 

The results are very encouraging. With 10 21 muon decays, the region Am^ > 5 x 
10~ 5 eV 2 is covered. For 5 x 10 21 muons, this region extends down to 10~ 5 cV 2 , in order 
words, the full range of values compatible with the LMA solar data is testable. 

If we consider that the value of Am^ is known and that it has a value of Am^ = 
10~ 4 cV 2 , one can constrain the values of the (5-phase within the range \5\ < 0.35 or 
\5\ > 2.8 for 10 21 muons and \5\ < 0.14 for 5 x 10 21 muon decays at the 90%C.L. We 
conclude that an exhaustive direct, model-independent exploration of the (5-phase, within 
the full range 10~ 5 < Am^ < 10~ 4 eV 2 requires an intensity of 5 x 10 21 muon decays of 
each sign. 

13 Conclusion 

We argued that the ideal neutrino detector at the neutrino factory should (1) have a mass 
in the range of 10 kton, (2) provide particle identification to tag the flavor of the incoming 
neutrino, (3) lepton charge measurement to select the incoming neutrino helicity, (4) good 
energy resolution to reconstruct the incoming neutrino energy, (5) possess a low energy 
threshold, in order to study neutrino events in the energy region where CP and T effects 
are the cleanest and most unambiguous (see Section 11) and (6) be isotropic to equally 
well reconstruct incoming neutrinos from different baselines (it might be more efficient to 
build various sources at different baselines, than various detectors). 

A detector with such qualities is most adapted to fully explore neutrino oscillations 
at the neutrino factory. In particular, measurement of the leading muons and electrons 
charge is the only way to fully simultaneously study CP and T violation effects. 

We think that a magnetized liquid argon imaging detector, scaled to 10 ktons, stands 
as the best choice of technique which holds the highest hope to match the above mentioned 
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Figure 14: Exclusion region at the 90%C.L. in the 5-phase vs Am^ plane. Two regions 
obtained with the Ay and Acp discriminants are shown. A 20% electron efficiency with a 
charge confusion probability of 0.1% has been assumed. The normalization is 10 21 (5 x 10 21 ) 
muon decays with energy En = 7.5 GeV and a baseline of L = 732 km. 



detector requirements. 

Such a detector could measure very precisely all the magnitudes of the elements of 
the mixing matrix and over-constrain them, because of its ability to reconstruct all final 
states including muons, electrons, tau-like and neutral current events (see Ref . JTT| ) . 

In addition, this detector could address the more challenging and most interesting 
goal of the neutrino factory, which is the search for effects related to the complex phase of 
the mixing matrix. This is because it could reconstruct the charge of both electrons and 
muons, and because it is perfectly adapted to reconstruct low energy events (typ. below 
15 GeV), which is the energy region where we expect these effects to be the cleanest and 
the most unambiguous. 

The choice of baseline and muon ring energy is in this context particularly critical. In 
view of the existence of massive devices like SuperK, MINOS or ICARUS, it is also worth 
considering if these detectors at their current baselines could be reused in the context of 
the neutrino factory. If new sites have to be found in order to satisfy the requirements of 
longer baselines, major new "investments" will be required. 

We find that the most favorable choice of neutrino energy E v and baseline L is in the 
region of the "first maximum" given by (L/E v ) max ~ 400 for |Ato 2 2 | = 3 x 10~ 3 eV 2 . 
This yields E™ ax ~ 2 GeV at 730 km, E™ ax ~ 8 GeV at 2900 km and E™ ax ~ 20 GeV 
at 7400 km for Am| 2 = 3 x 10~ 3 eV 2 . 

We showed that the discriminants of the <5-phase all scale with L/E v . This property 
leaves flexibility in the choice of the actual neutrino energy and the baseline, since only 
their ratio L/ E v is determinant. Because of the rising neutrino cross-section with energy, 
it will be more favorable to go to higher energies if the neutrino fluence is constant. 
Keeping the L/E v ratio constant, this implies an optimization at longer baselines L. 
One will hence gain with the baseline L until we reach L max « 5000 km beyond which 
matter effects will effectively reduce the dependence of the oscillation probabilities on the 
(5-phasc. 

As far as the neutrino energy is concerned, it is clear that the average neutrino energy 
E v scales linearly with the muon beam energy E^ . A non- negligible aspect of the neutrino 
factory is the need to accelerate quickly the muons to the desired energy, and so, it is 
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expected that higher energies will be more demanding that lower ones. Eventually, cost 
arguments could determine the muon energy. It could therefore be that lower energy, more 
intense neutrino factories could be more advantageous than higher energy, less intense 
ones. 

From the above arguments, we therefore conclude that a very intense neutrino factory, 
capable of providing more than 10 21 useful 7.5 GeV muon decays, directed towards a 
distance of 730 km, coupled with a magnetized liquid argon imaging detector of about 
10 kton would provide the ultimate, most comprehensive setup to study CP and T 
violation effects in neutrino flavor oscillation. 
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